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Abstract 

The recent work of Delabaere and Trinh (2000 J. Phys. A 33 8771) discovered 



> 

cn 

O ■ o 

Tjlj- . the existence of "PT-symmetry breaking, complex energy, L solutions for the 



one dimensional Hamiltonian, P^ + iX^ + iaX, in the asymptotic limit, a — > 
— oo. Their asymptotic analysis produced questionable results for moderate 



o 

r~| ■ values of a. We can easily confirm the existence of "PT-symmetry breaking 



solutions, by explicitly computing the low lying states, for \a\ < 0(10). Our 
analysis makes use of the Multiscale Reference Function (MRF) approach, 
developed by Tymczak et al (1998 Phys. Rev. Lett. 80 3678; 1998 Phys. 
Rev. A 58, 2708). The MRF results can be validated by comparing them 
with the converging eigenenergy bounds generated through the Eigenvalue 
Moment Method, as recently argued by Handy (2001a, b). Given the reliability 
of the MRF analysis, its fast numerical implementation, high accuracy, and 
theoretical simplicity, the present formalism defines an effective and efficient 
procedure for analyzing many related problems that have appeared in the 
recent literature. 

PACS numbers: 02.30.Hq, 03.65.-w, 03.65. Ge 
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I. INTRODUCTION 



There has been much interest, recently, in understanding the symmetry breaking mech- 
anism for invariant Hamiltonians of the type + Y,j=oCj{iXy . The mathematical 
interest in these systems originated from a conjecture by D. Bessis, and WKB analysis confir- 
mation by Bender and Boettcher (1998), that the class of potentials of the form V{x) — (iX)" 
only allow for 'PT-invariant solutions, and thus can only have real discrete spectra. The 
recent literature testifies to the great interest in these problems, as can be found in the cited 
references by Bender et al (2001), and Mezincescu (2000, 2001). 

An important work establishing that PT-symmetry breaking systems do exist was the 
recent study by Delabaere and Trinh (2000) which used asymptotic methods to analyze 
the Hamiltonian Ha = + iX^ + iaX. Their analysis showed the existence of symmetry 
breaking solutions for large a values; however, for moderate values, their results were near the 
limits of their analytical validity. Despite this, they made certain questionable predictions 
for moderate a values. 

Our immediate objective is to check the validity of the Delabaere and Trinh study, by 
explicitly computing the low-lying complex (and real) eigenenergies for moderate a values. 
We do this in two ways. The first makes use of the very efficient Multiscale Reference 
Function (MRF) formalism of Tymczak et al (1998a,b). The results of this eigenenergy esti- 
mation analysis are then confirmed through application of a recently developed eigenenergy 
bounding theory proposed by Handy (2001a,b), which can generate converging bounds to 
the complex eigenenergies. This bounding approach, referred to as the Eigenvalue Moment 
Method (EMM), is exact, although numerically slower in its implementation than the MRF 
procedure. 

As a footnote to the above, we emphasize that the EMM theory generates an infinite 
hierarchy of closed, finite dimensional, algebraic, eigenenergy constraints. These are then 
solved, numerically. In this regard, the EMM procedure is very different from other nu- 
merical schemes, such as numerical integration, which are intrinsically of an approximating 



2 



nature, and cannot provide any fundamental theoretical insight into the underlying physical 
processes. The algebraic constraints generated through the EMM formahsm can provide 
such insight; although such analysis have not been attempted, heretofore. However, for the 
immediate purposes of this study, we solely defer to EMM in order to check the validity of 
the MRF results. 

In this work, we provide the essentials of the MRF theoretical structure, as applied to 
the Ha Hamiltonian. The EMM theory is not discussed. Only the numerical bounds are 
quoted in the Tables. This work vahdates the relevancy of MRF theory in the computation 
of complex eigenenergies, for the class of problems referenced above. 

Both the MRF and EMM methods are dependent on a moments' representation for the 
given system. This in turn is readily realizable for any (multidimensional) rational fraction 
potential. 

Any moment based analysis is inherently multiscale in nature. That is, as the number 
of moments used increases, one is probing the system at successively smaller scales. Consis- 
tent with this, the MRF basis representation, particularly within configuration space, has 
important ties with (complex) turning point quantization (Handy et al (2000)), and wavelet 
analysis (Handy and Brooks (2001)). 

We outline the basic MRF theory and its implementation, in the next section. The last 
section contains a detailed enumeration, and illustration, of the MRF results, which make 
precise the qualitative spectral structure conjectured by Delabaere and Trinh. 
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II. THE MRF REPRESENTATION 



A. The Moment Equation 

The starting point for the MRF analysis is the transformation of the Schrodinger equation 
into the Fourier representation, assuming that one is working with the physical, L^, solutions. 
Thus, for the configuration space Schrodinger equation studied by Delabaere and Trinh 
(2000) 

-dl'^{x) + {ix^ + iax)^{x)^E^{x), (1) 
its Fourier transform counterpart is (i.e. x — > idk, ik) 

k''^{k) + {dl-adk)'l/{k) = E^{k), (2) 

where 

^{k) = dx e-'^''^{x). (3) 

It is important to note that, for the physical solutions, a simple apphcation of WKB 
analysis (Bender and Boettcher (1998)) tells us that the asymptotic behavior of the configu- 
ration space representation yields an entire Fourier transform. Because of this, the A;-power 
series expansion is absolutely convergent, and defined in terms of the power moments: 

where 

/■ + OC 

//p = / dx x^'^{x), (5) 



define the Hamburger power moments. 

We can generate the recursion relation for the /Xp's from the standard power series ex- 
pansion methods for linear differential equations (Bender and Orszag (1978)). Alternatively, 
we can apply dx x^ to both sides of Eq.(l), combined with integration by parts, and 
obtain the necessary Moment Equation: 



Atp+3 = -o^^^p+l - iEiip - ip{p - l)Atp-2, (6) 

for p > 0. This corresponds to a homogeneous, hnear, finite difference equation, of effective 
order 1 + = 3, since specification of the independent moments {fJ-o, fJ-i, ^J^2}, plus the 
(complex) energy parameter, E, generates all of the remaining moments. The independent 
moments are referred to as the missing moments. 

The linear dependence of the moments, on the missing moments, can be expressed 
through the relation 

lip = Y;^Mp^e{E)i^e, (7) 

e=o 

where the energy dependent coefficients satisfy the moment equation, with respect to the 
p-index, as well as the initial conditions Mi^^£^ — Sj^^^^^^, for < £i,i2 < ms- 

B. Defining an Analytic Basis in the Fourier Space 

The physical solutions in the Fourier representation must also be L^. One would like to 
find an appropriate basis into which to transform the Fourier power series expansion: 

4=ti-^kr^ = 4=taMk). (8) 

V27r „=o P- V ^TT 



p=0 -f^- V j=o 

The easiest choice, leading to a rapid, analytic generation of the aj coefficients, is to take 

Bj{k) = {-ikyR{k), (9) 

where R is some arbitrary "reference" function yielding a complete (if not orthogonal) basis. 
Thus, if is analytic, then one can generate the a^'s by expanding In particular, we 

can take R{k) = e"^*^^, where /5 > 0, and otherwise arbitrary, yielding 



oo 



A{k) ^j:a,{-zky = e^'' EHkf^- (10) 

It is then clear that the a/s will become linear in the missing moments, {//^lO < i < rUs}. 
Specifically, 



j{E; ^0, 1^1, 1^2) = 



p+2q=j 



q\p\ ' 



or 



rns=2 



aj(E; iiQ, III, 1^2) ^ J2 { 



e=0 p+2q=j 1-P- 



(11) 



(12) 



for j > 0. Clearly, the maximum (Hamburger) moment order generated, Pmax, determines 
the maximum order of a/s generated, < j < Pmax- 

C. The MRF Quantization Prescription 



It has been argued by Tymczak et al (1998a,b) that the convergent zeroes of the coefficient 
functions 



converge to the exact discrete state energies: 

lim E^^^ = ^P%sicai_ 



(13) 



(14) 



The above root equation must be adapted to the, l+rUg, linear, missing moment structure 
of the Hamiltonian in question. Thus, to any expansion order J, we impose that the last 
1 + 777.5, Oj-coefficients be zero (i.e. aj-i = 0, for < £ < rUg). This results in a 1 + ms = 3 
dimensional, determinantal equation for the energy: 



where 



Aj{E) = Det 



Aj{E) = 0, 



/A<^(E, /3), mS(E, /3), -Mg(i?, /?) \ 

{E, (3) , m['} {E, (3) ,M^i}{E, P) 
\Mi',^(E,(3),Mi'}(E,(3),Mi',^{E,(3) J 



(15) 



(16) 



and 



M\i;eAE,(3)^ Y 

p+2q=J-h 



'-P)m,,e,{E) 
q\p\ 



(17) 



for < £1, £2 < = 2. In the Tables, the Pmax parameter corresponds to J = P,^ 



III. NUMERICAL IMPLEMENTATION OF MRF 



The MRF analysis is implemented for /? = .5. This value of the arbitrary parameter 
generates the fastest converging results. 

In Figure I we plot the PT-symmetry breaking solutions {Im{E) ^ 0), and the VT- 
symmetry invariant solutions {Im{E) — 0), for moderate a values at the limits of Delabaere 
and Trinh's asymptotic analysis. There are four branches depicted. In Tables I- IV we specify 
some of the points plotted, and in addition, compare the MRF eigenenergy estimates with the 
EMM bounds. The MRF results in these tables generally correspond to J < 50. The EMM 
results were calculated to the same moment order, Pmax < 50, for the complex energies. In 
those cases where the MRF method predicts a purely real eigenenergy, we used the faster 
EMM formahsm corresponding to exphcitly assuming that the underlying solutions are VT 
invariant (Handy (2001a)). For this case, the underlying problem is of Stieltjes character, 
and a Stieltjes moment order of Pj^^^ = 30, corresponds to a Hamburger moment order of 60 
(i.e. Pjnax — ^Y^)- For this reason, in some cases the MRF estimates (computed at J < 50) 
he outside of the EMM bounds. However, for some of the EMM bounds, numerical instabihty 
concerns required that P^^ax < 30 (these were computed on a standard, double precision, 
IBM platform). Such cases are relatively few in number, and are explicitly identified. By 
working at a larger precision order, the corresponding bounds can be improved. 

It will also be noted that for a = 0, we quote the EMM bounds generated at higher 
Stieltjes moment order {Pj^^^ < 60), given in the work by Handy (2001a). These bounds 
are very consistent with the corresponding entries in Tables V-VI. 

In Tables V and VI we quote the MRF results, for J = 100 (using CRAY-double preci- 
sion). Only the stable digits are given (that is. Tables V and VI correspond to the (empiri- 
cally determined) stable digits within the MRF generated sequence, for < J < 100). 

The results are consistent, that is, the MRF results in Tables I-IV, for the most part, 
lie within the bounds. This is definitely the case for Tables V - VI. Note that in the Tables 
we quote the imaginary parts of the energy as ±, this is because we cannot tell which VT- 



7 



invariant branch continues into the T'T-breaking branch. For VT invariant Hamiltonians, 
complex energies come in conjugate pairs. Thus, if £J is a solution, so too is E* . 

In Figures II - III we narrow in on the smaller of the two critical a values, acn ^) as 
shown in Fig. I. At these critical points, the energy goes from being real {VT - invariant 
solutions) to complex [VT - breaking solutions) . They are: 

a^^ = -2.611809356, (18) 
corresponding to E'en = 1.28277353562; and 

a^^ = -5.375879629, (19) 
corresponding to £^^2 = 4.181388093. 
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TABLES 



TABLE I. £^i-Branch for the Discrete States of + iX^ + iaX 

a E^MRF) EMM Bounds : E^^^ < Er < E^ h Ef^ < Ej < E^ 

-5.0 ( 1.3433409, ± 2.9073602) 1.343311^ < Er < 1.343354\ 2.9073^ < Ej < 2.9075^ 

-4.5 ( 1.2992519, ± 2.3124924) 1.299242 < Er < 1.299252, 2.3124 < Ej < 2.3126 

-4.0 ( 1.2486637, ± 1.7617076) 1.248637 < Er < 1.248666, 1.761688 < Ej < 1.761742 

-3.5 ( 1.2124399, ± 1.2609114) 1.212421 < Er < 1.212448, 1.26088 < Ej < 1.26094 

-3.0 ( 1.2258438, ± 0.7600296) 1.225837 < Er < 1.225864, .76000 < Ei < .76004 

-2.5 ( .9280136, ) 0.92799980^ <E< 0.92800101^ 

-2.0 ( .6209137, ) 0.62091347 <E< 0.62091386 

-1.5 ( .5964936, ) 0.59649326 <E < 0.59649351 

-1.0 ( .6999615, ) 0.69995977 <E< 0.69995978 

-.5 ( .8926699, ) 0.89266827 <E< 0.89266849 

.0 ( 1.1562673, ) 1.15626695 <E < 1.15626718 

1.156267071988113243 <E< 1.15626707198811335=^ 

.5 ( 1.4798519, ) 1.47985179 <E< 1.47985206 

1.0 ( 1.8561128, ) 1.85611065 <E < 1.85611108 

1.5 ( 2.2797563, ) 2.27975185 <E< 2.27975232 

2.0 ( 2.7467434, ) 2.74673952 <E< 2.74674023 

2.5 ( 3.2538767, ) 3.25387596 <E < 3.25387723 

3.0 ( 3.7985559, ) 3.79855387 <E< 3.79855395 

3.5 ( 4.3786140, ) 4.37859645 <E< 4.37859736 

4.0 ( 4.9921974, ) 4.99215436 <E < 4.99215504 

4.5 ( 5.6376822, ) 5.63763149 <E< 5.63763200 

5.0 ( 6.3136428, ) 6.31359739 <E< 6.31360665 



^EMM Analysis of Handy (2001b) 

9 



^EMM {Pj^l = 30) Analysis of Handy (2001a) 
^EMM (P^al = 60) Analysis of Handy (2001a) 
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TABLE II. £;2-Branch for the Discrete States of + iX^ + iaX 



a 



EMM Bounds : E^^^ < Er < E^^ k Ef^ < Ej < e\ 



-5.0 
-4.5 
-4.0 
-3.5 
-3.0 
-2.5 
-2.0 
-1.5 
-1.0 
-.5 
.0 

.5 
1.0 
1.5 
2.0 
2.5 
3.0 
3.5 
4.0 
4.5 
5.0 



( 1.3433409, ± 2.9073602) 
( 1.2992519, ± 2.3124924) 
( 1.2486637, ± 1.7617076) 
( 1.2124399, ± 1.2609114) 
( 1.2258438, ± 0.7600296) 
1.6859358, 
2.2922626, 
2.7425268, 
3.1797220, 
3.6320373, 
4.1091279, 



4.6147402 
5.1501688 
5.7154576 
6.3100192 
6.9332376 
7.5850094 
8.2656580 
8.9745543 
9.7077326 



( 10.4575130, ) 



1.343311^ <Eb.< 1.343354\ 2.9073^ < Ei < 2.9075^ 

1.299243 <Er< 1.299252, 2.3124 < Ei < 2.3126 
1.248637 < < 1.248666, 1.761688 < Ei < 1.761742 
1.212421 <Er< 1.212448, 1.26088 < Ej < 1.26094 
1.225837 <Er< 1.225864, .76000 < Ei < .76004 
1.68597765^ <E< 1.68598087^ 
2.29229055 <E < 2.29229333 
2.74252667 <E < 2.74253034 
3.17971312 <E< 3.17971750 
3.63207237 <E< 3.63207767 
4.10922704 <E < 4.10923558 
4.109228752806^ <E < 4.109228752812=^ 
4.61483391 <E < 4.61484633 
5.15016059 <E < 5.15017640 
5.71538438 <E< 5.71541649 
6.31020527 <E < 6.31025282 
6.93405453 <E < 6.93412954 
7.58627841 <E< 7.58638028 
8.26613065 <E < 8.26633647 
8.97272640^ <E < 8.97332836^ 
9.70210565^ <E< 9.70789436^ 
10.45227656^ <E< 10.48006875^ 



^EMM Analysis of Handy (2001b) 

^EMM {P^}^ = 30) Analysis of Handy (2001a) 
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=^EMM (P^^2. = 50) Analysis of Handy (2001a) 
^EMM (P^al < 30) Analysis of Handy (2001a) 
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TABLE III. E^a-Branch for the Discrete States of + iX^ + iaX 



Cm 


-^3 I"/ 


EMMRmmHs ■ E^^^ < Eo < E^'^^ Rr R*-"^^ < Et < E^^^ 


7 ^ 




A Ai~)fi.fif^9^ ^ P,> ^ A A^ 149(^1 9fizL'^ ^ Pr ^ 9700 




l' 4 ^19n^/lQ + 9 'ifi'^47'i1'l 




~\J.O 


/' 4 9n^fi^8n + 1 Qfl^Qfi'iS'l 
I ^.zuouooUj in ±.e7UOc/uooj 


4 904^90 ^ ?^?o <r 4 90(^^(^0 1 009*^ <^ Kt <^ ^ 0074 


-u.u 




A ^'^'\R(^A <r Rr> <r A 1'^4^'^4 1 9RA <^ -\-F,t ^ ^ ORf^ 








-o.yj 




^ 4^1 ^R741 2 ^ ^ Q 4Q1 ^Qi^402 




( ^ ^9^944fi n "1 


q QOQOQl^QC; ^ F ^ ^ ^9^91^01 S 


A. 




■^0k7R0QQ ^ F ^ ^ '^0k77'^44 


~o.o 


I 0.01 (Ue70±j U J 


^ R77f{RdR^ <^ R <^ 877701 7R 

O.Ol 1 UO^OO \ J-J \ O.Ol 1 lUllO 


~o.u 




4 ^^^4997^ ^ ^ 4 ^^^44^1^4 






4 k99Q4R4Q ^ F ^ 4 899Q8497 
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^ ^^l^'=i'^^7 ^ R ^ ^ ^^140^11 

O.OOlOOOOf \. J-J \ 0.001'4:U01± 


-1.0 




O.OOlOyOlO \ J-/ <v.. O.OOlUUOlU 


1 


( f\ 4n'ifi7f^9 n "1 


R 40^^89^10 ^ F ^ R 40^79RQ4 






R Q71 ^^Qi^l <r R <^ f\ Q71 "^97^^ 


n 
.u 




7 "^^91 "^QOl <r R <^ 7 '^R949^'=»'=» 






7 "^8997^8^^403 ^ R ^ 7 '^R997'^8'^'^1 ^ 


.0 


(" 8 177014^ n ^ 


R 1 7R87901 <^ R <^ R ^ 7790R44 


J- .u 


^' 8 81 f>9941 n "1 


8 81 "^iRO^Rl <^ P ^ 8 81 R44^7'=i 


±.0 


I q 480901 1 n "1 


Q 47Q1 ^"^04 <' R <' Q 480449^0 


2.0 


( 10.1686764, ) 


10.16683597 < E < 10.16907546 


2.5 


( 10.8807993, ) 


10.87896180 < E < 10.88253390 


3.0 


( 11.6159802, ) 


11.61535000 <E < 11.62075000 


3.5 


( 12.3755952, ) 


12.37090000^ <E< 12.38980000^ 


4.0 


( 13.1638132, ) 


13.16160000^ <E< 13.20480000^ 
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4.5 
5.0 



( 13.9842043, ) 
( 14.8298186, ) 



13.86000000^ <E< 14.04000000^ 
14.72000000^ <E < 15.08000000^ 



^EMM Analysis of Handy (2001b) 
2EMM {Pl^^fl = 30) Analysis of Handy (2001a) 
3EMM (P^al = 50) Analysis of Handy (2001a) 
^EMM (P^l < 30) Analysis of Handy (2001a) 
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TABLE IV. £;4-Branch for the Discrete States of + iX^ + iaX 



a Ef^^^^ (a) EMM Bounds : E^^^ <Er< E^^^ k Ef^ < Ej < Ef^ 

-7.5 ( 4.4103093, ± 3.2664292) 4.408662^ <Eb.< 4.411426\ 3.2645 < Ei < 3.2700 

-7.0 ( 4.3120349, ± 2.5634751) 4.310530 <Er< 4.313107, 2.5618 < Ei < 2.5674 

-6.5 ( 4.2056380, ± 1.9039658) 4.204520 <Er< 4.206360, 1.9026 < Ei < 1.9074 

-6.0 ( 4.1343003, ± 1.2850914) 4.133864 <Er< 4.134534, 1.284 < ±Ei < 1.286 

-5.5 ( 4.1575135, ± 0.5313291) 4.155648 < Er < 4.158235, .5300 < Ej < .5370 

-5.0 ( 5.1678291, ) 5.16784214^ <E< 5.16798149^ 

-4.5 ( 5.8041737, ) 5.80427900 <E < 5.80439660 

-4.0 ( 6.3796826, ) 6.37969997 <E < 6.37985480 

-3.5 ( 6.9490904, ) 6.94880880 <E< 6.94915872 

-3.0 ( 7.5252398, ) 7.52489539 <E< 7.52536568 

-2.5 ( 8.1130836, ) 8.11302413 <E < 8.11368768 

-2.0 ( 8.7159047, ) 8.71631657 <E< 8.71698842 

-1.5 ( 9.3364444, ) 9.33658547 <E< 9.33774669 

-1.0 ( 9.9763275, ) 9.97552896 <E < 9.97718784 

-.5 ( 10.6352870, ) 10.63343056 <E< 10.63660576 

.0 ( 11.3120046, ) 11.31165120 <E< 11.31588480 

11.314421818^ <E < 11.314421824^ 

.5 ( 12.0072541, ) 12.01004000 <E< 12.01676000 

1.0 ( 12.7264921, ) 12.72473600 <E < 12.74288000 

1.5 ( 13.4761039, ) 13.47000000 <E < 13.51500000 

2.0 ( 14.2487931, ) 14.21000000 <E< 14.33000000 

2.5 ( 15.0088412, ) 14.96800000 <E< 15.22000000 

3.0 ( 15.8347839, ) 15.71250000 <E < 16.12500000 

3.5 ( 16.4020563, ) 16.40000000 <E< 17.00000000 
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^EMM Analysis of Handy (2001b) 

^EMM {P^l = 30) Analysis of Handy (2001a) 

^EMM {P^^^ = 50) Analysis of Handy (2001a) 
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TABLE V. MRF Complex Eigenenergy Estimates for P^ax = 100 



a 


Ei&z, E2 




-7.5 




(4.4102443440, ± 3.2664976653) 


-7.0 




(4.3119949171, ± 2.5635584867) 


-6.5 




(4.2055983584, ± 1.9040248046) 


-6.0 




(4.1342519473, ± 1.2851227083) 


-5.5 




(4.1574750374, ± 0.5313567960) 


-5.0 


(1.3433431987, ± 2.9073906160) 




-4.5 


(1.2992423296,± 2.3125154783) 




-4.0 


(1.2486567335,± 1.7617193016) 




-3.5 


(1.2124367296,± 1.2609099725) 




-3.0 


(1.2258475767,± 0.7600224714) 
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TABLE VI. MRF Real Eigenenergy Estimates for P^ax = 100 



a 



E2 



E3 



E4 



-5.0 
-4.5 
-4.0 
-3.5 
-3.0 
-2.5 
-2.0 
-1.5 
-1.0 
-0.5 



.5 
1.0 
1.5 
2.0 
2.5 
3.0 
3.5 
4.0 
4.5 
5.0 



.928000342158 
.6209135740648 
.59649338409590 
.6999599207986 
.8926684335546 
1.1562670719881 
1.47985186079689 
1.8561107660566 
2.2797520475930 
2.7467399808560 
3.2538769263689 
3.798554700716 
4.37859694562365 
4.9921540825786 
5.637630445616 
6.313632040368 



1.685979342731 
2.2922925019 
2.7425293939422 
3.1797157763716 
3.632074461334 
4.109228752809 
4.6148387273616 
5.150168955614 
5.71540870715 
6.31023836105 
6.934096040540 
7.586310988692 
8.2661728690 
8.9729684347 
9.7060008600 
1.046459957615 



3.4313832016 
3.3232426191 
3.5087656073 
3.8776930052 
4.3334398364 
4.822975394502 
5.3313834010 
5.8576220019 
6.4036451049 
6.971403910342 
7.5622738549 
8.1770825316 
8.8162451717 
9.47988920956 
10.1679478499 
10.880225999 
11.616445659 
12.37627693 
13.15935916 
13.96531533 
14.7937618 



5.1678886857 
5.8043392228 
6.3798052063 
6.9490523010 
7.5251919556 
8.113377225 
8.716619963 
9.337028976 
9.97613238 
10.63501626 
11.3144218 
12.01482435 
12.73649738 
13.4795632 
14.2440319 
15.02983115 
15.8368284 
16.6648481 
17.513684 
18.383109 
19.2728823 
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IV. CONCLUSION 



We have confirmed the asymptotic analysis prediction of Delabaerre and Trinh (DT) 
on the existence of symmetry breaking solutions for the Hamiltonian. Our methods 
enable the precise analysis of the complex-real spectra, particularly for moderate a values 
at the hmits of their (DT) asymptotic vahdity. The results of both an eigenenergy estimation 
method (MRF) and an eigenenergy bounding method (EMM) were presented. The algebraic 
simplicity, and ease of computational implementability, of the MRF method recommend it 
highly for application to similar problems. Through the use of readily available algebraic 
programming software, the MRF approach can be extended to arbitrary precision (indeed, 
Tymczak et al (1998b) were able to generate the quartic anharmonic oscillator ground state 
energy to more than 250 decimal places), making it a very powerful tool in these types of 
investigations. 
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